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THEORY  OF  OCEAN-ARRAY  INTERACTION, 

PHASE  I  REPORT 

1.  INTRODUCTION 

A  linear  model  describing  the  transverse  dynamics  of  a  long  flexible  cylinder  in 
axial  flow  is  presented  in  references  1  and  2.  Experimental  verification  of  the  theory 
is  discussed  in  reference  3.  That  work  describes  the  transverse  response  of  the 
cylinder  to  a  forced  vibration  at  the  upstream  end  of  the  cylinder. 

The  analytical  procedure  of  the  previous  work  involved  solving  a  homogeneous 
momentum  equation  in  which  the  forced  vibration  entered  through  the  upstream 
boundary  condition.  The  solution  modeled  a  cable/towed  array  system  excited  by 
the  random  crosstrack  meandering  of  the  towing  vessel.  In  this  study  we  are  con¬ 
cerned  with  the  transverse  dynamics  of  a  long  flexible  cylinder  in  “nearly”  axial 
flow.  In  the  present  case,  a  space-time  field  of  fluid  flow  orthogonal  to  the  axial 
flow  is  superimposed  on  the  basic  axial  flow.  This  crosscurrent  is  a  perturbing  flow 
that  always  has  a  magnitude  much  smaller  than  the  axial  flow.  The  intent  of  this 
study  was  to  model  the  transverse  dynamics  of  a  cable/array  system  driven  by 
crosscurrents. 

Analytically,  the  procedure  in  the  present  work  involved  solving  the  same 
previously  discussed  momentum  equation,  only  now  modified  to  be  non- 
homogeneous  with  the  oceanic  excitation  entering  through  the  nonhomogeneous 
term. 


The  nonhomogeneous  governing  equation  is  derived  and  discussed  in  section  2. 
Approximate  solutions  relevant  to  low-  and  high-frequency  domains  of  interest  are 
discussed  in  section  3.  In  section  4,  a  complete  solution  is  developed  and  a 
numerical  technique  is  applied  to  the  problem  in  section  5.  The  work  is  summarized 
in  section  6.  Excitation  of  the  system  by  a  plane  wave  is  discussed  in  each  section. 


2.  GOVERNING  EQUATION 
2.1  MOMENTUM  EQUATION 

The  basic  partial  differential  equation  (PDE)  applied  in  this  study  was  derived 
and  discussed  in  several  papers.  The  equation  describes  the  dynamics  of  a 
physical  configuration  (see  figure  1).  The  configuration  consists  of  a  flexible 
cylindrical  body  of  circular  cross-section  that  is  immersed  in  an  incompressible 
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fluid  of  uniform  density.  The  fluid  is  flowing*  with  uniform  velocity  (U)  parallel  to 
the  x  axis. ^  The  equation  is  expressed  as 

(M  +  ■£-§  +  mo2^-|  +  2MU^£  -  JL02(l  .  ,)£ 

-°  • 

where 

y(x,t)  a  crosstrack  position  of  the  cylinder  measured  at  x, 
x  s  alongtrack  independent  space  variable  increasing  downstream, 

t  s  time, 

M  3  effective  mass  of  fluid  “pushed”  by  the  cylinder  per  unit 
length  of  cylinder, 

m  3  mass  of  the  cylinder  per  unit  length, 

U  3  free  stream  velocity, 

L  3  cylinder  length, 

dc  3  cylinder  diameter, 

ct  3  tangential  (longitudinal)  drag  coefficient,  and 
cn  3  coefficient  of  the  linearized  normal  drag  expression. 


•Or,  equivalently,  being  towed  with  uniform  velocity  (U). 
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y*  AXIS 


The  effetfof  a  current^  component,  V,  in  the  y  direction  is  to  modify  the  original 
flow  from  Ui  to  ut  +  Vj,  where? and  ? are  unit  vectors  in  the  coordinate  system  of 
figure  1.  This  flow  modification  is  important  when  examining  the  hydrodynamic 
terms.  It  is  important  to  recognize  that  because  hydrodynamic  forces  are 
nonlinearly  related  to  fluid  velocities,  we  cannot  calculate  the  Cartesian  components 
of  the  force  directly  from  the  Cartesian  coordinates  of  the  fluid  velocity.  We  must 
first  find  the  total  fluid  velocity  and  direction  and  then  the  resultant  vector  force. 
Figure  2  shows  the  geometric  relationships.  Note  that  the  figure  ignores  the  small 
angle  assumptions  inherent  in  the  derivation  of  equation  (1)  in  order  to  improve  the 
geometric  visualization. 

The  cylinder  is  towed  in  the  negative  x  (Lection  creating,  in  the  translating 
coordinate  systems  of  figures  1  and  2,  a  flow  component,  U,  in  the  positive  x 
direction.  The  crossflow,  V,  is  in  the  positive  y  direction  which  creates  a  local  total 
fluid  velocity  having  polar  coordinates  U  [1  +  (V/U)2]1/2  and  $.  This  is  a  space- 
time  function  because 


a  A, 

V  -  V(x,t)  . 

The  normal  hydrodynamic  drag  on  the  cylinder  is  well  known  to  be 

m=”2[i  4T]fta)  • 

where  er  is  the  angle  between  the  fluid  velocity  vector  and  the  tangent  to  the  cylinder 
at  a  point.  The  form  of  f(o)  used  in  reference  1  and  later  experimentally  supported 
in  reference  6  may  be  written  as 
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!c  sin  a  for  a  <.  3°  . 

11  2 

c*  sin  a  for  a  >  3°  . 
n 

Typical  numerical  values  of  these  coefficients^  are  cn  =  1.2  and  c^  =  0.07.  A 
basic  assumption  of  equation  (1)  is  that  cr<  3*.  Here  we  are  further  going  to  require 
that 


A 


V/U  «  1  . 


Thus,  the  normal  hydrodynamic  drag  becomes 

l/2pd  c  U2a  for  a  <  3°  .  (2) 

c  n 

Lighthill7  showed  that  the  fluid  velocity  relative  to  the  cylinder  and  normal  to  the 
axis  of  the  cylinder  is 


hi 

3t 


+  U 


hL 

dx 


9 


A  AA 

where  y'  and  x  are  coordinates  determined  b^the  total  fluid  flow,  Ui  +  Vj,  and  are 
defined  in  figure  2.  We  again  assumed  that  V/U  «  1 .  Following  Paidoussis,4  we 
make  the  small  angle  approximation* 


a  st 


+  U 


U 


(3) 


and  introduce  the  variable 


The  transformation  from  the  primed  to  unprimed  coordinate  system  is  shown  in 
figure  2.  x0,  y0  is  an  arbitrary  position  and  can  be  taken  to  be  zero.  Following  the 

rules  of  partial  derivatives  we  may  write  the  following  relations: 

!£-  -  -<x  +  y  tan  40  f|  -  sin  4>  +  f* 
and 

4^-  -  -(x  +  y  sin  (J>)  |£  -  sin  d>  ||  +  cos  <|>  |^-  . 


*Note  that  a  is  expressible  as  the  arcsine  of  the  normal  and  axial  components  of 
flow  velocity. 


4 


The  small  angle  assumption  allows  the  following  relations:  cos  =  1  and  sin  $  =  <j>  = 
tan  $  and  dx/dt  2  0.  Thus  the  partial  derivatives  become 


and 


!£-■  &  -  *  +  [(*  +  If] 


(4a) 


(4b) 


Because  the  transformations  are  “local”  (i.e.,  they  are  applied  at  x  =  y  =  0),  we 
may  substitute  the  expressions 


hL 

3x' 


A 

1 

3x  "  U 


and 


at  St 


into  equation  (3).  This  makes  the  normal  hydrodynamic  term  equal  to 


KCnU  (lf  +  U|j)  ’ 


A 

uv 


We  may  modify  this  equation  by  eliminating4^  c  using  the  form 

p  =*  4M/ird*  , 

yielding 


2cn  MU  / 3y  ,.3y  \  2cn  MU  * 

tt  dc  \3t  ^x )  v  dc 


(5) 


The  first  term  above  already  appears  in  equation  (1).  The  second  term  is  the 
nonhomogeneous  addition  and  becomes  the  right-hand  side  of  our  present  problem. 

Note  that  the  presence  of  V  in  the  problem  does  not  alter  the  inertial  terms  of 
equation  (1)  and,  assuming  V/U«  1,  does  not  alter  the  tension  term  either. 


The  following  nondimensional  variables  are  introduced  into  equation  (1)  modified 
by  equation  (5): 

A 


X 


tu 

L 


0 


M 

M  +  m 


E  -  j-  ,  and  V  -  1  ,  (6) 


c 


where  A  is  an  arbitrary  length  which  can  take  on  any  value  because  equation  (1)  is 
linear  in  y.  The  result  is 


5 


where 


a  «  8(1  -  2cte)  , 
b  -  28cte  , 

c  ■  2(c  +  c  /7r)eg  ,  and 

t  n 

d  *  (2/tt)8c  e  . 

n 

The  dependence  of  the  nonhomogeneous  term  of  equation  (7)  on  cn  is  intuitively 
clear,  however,  the  dependence  on  z  is  not  as  obvious.  The  role  of  z  in  that  term  is  to 
show  the  relative  importance  of  the  nonhomogeneous  hydrodynamic  forcing 
function  relative  to  the  inertial  forces,  i.e.,  the  first  two  terms  of  equation  (7).  Large 
values  of  z  reduce  the  inertial  effects,  ft  for  a  cylinder  at  relatively  small  frequencies 
is  a  constant  (i.e.,  1/2). 

2.2  BOUNDARY  CONDITIONS 

The  appropriate  boundary  conditions  have  been  discussed  in  references  1,  4,  5, 
and  8.  Two  boundary  conditions  are  required  for  the  second-order  description  of 
equation  (7).  In  general,  the  upstream  (4  =  0)  boundary  condition  is  an  arbitrary 
time  function,  that  is, 

n<5  -  o , t)  -  g(x)  . 

For  a  boundary  condition  on  the  downstream  end  of  the  cable  (4  =  1),  we  en¬ 
counter  a  fundamental  problem  that  is  a  consequence  of  dropping  the  flexural 
rigidity  term  in  the  Paidoussis  momentum  equation.  Because  the  physical  constraint 
at  the  free  end  involves  bending  mr  ments,  the  second-order  system  (which  allows  no 
bending  moments)  is  unable  to  characterize  the  free  end  physical  condition. 
Analytically,  this  problem  results  from  changing  from  a  fourth-order  to  a  second- 
order  description. 

Historically,  this  type  of  problem  (D.  Bernoulli’s  “hanging  chain,”  1732)  has 
been  solved  by  invoking  a  boundedness  condition,  that  is, 

jn(C,T)|  <  oo  for  0  <.  £  £  1  . 

It  will  be  useful  to  think  of  the  nondimensional  displacement,  rj,  as  consisting  of 
two  additive  terms,  that  is, 

n  -  Oj  +  n2  . 

Thus,  with  the  linear  operator  D  defined  by  equation  (7)  we  may  write 
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with 

n(5.T)  -  g(t) 

for  £  -  0 

and 

|n(£,T) |  <  « 

for 

0  <  £  <  1  , 

or,  equivalently,  we  have  the  system  of  equations 

Dp^C.t)  -  0 

for 

0  <  £  <  1  . 

n^C.T)  -  g(T) 

for 

£  -  0  . 

In^.T)!  <  » 

for 

0  <  £  <  1  . 

and 

Dri2(£,T)  -  -dV(5,T) 

for 

0  <  £  <  1  . 

ti2(£,t)  -  0 

for 

£  -  0  . 

|n2(£,T) |  <  « 

for 

0  <  £  £  1  . 

(7b) 


(7c) 


Thus  r/i  and  rjj  are  solutions  to  two  types  of  problems:  (1)  ni  is  the  solution  to  the 
homogeneous  momentum  equation  and  is  the  response  to  the  forced  vibration 
problem  discussed  in  references  1  and  2,  and  (2)  >72  is  the  solution  to  the 
nonhomogeneous  equation  and  is  the  response  to  the  crossflow.  Note  that  in  general 
172  will  equal  the  sum  of  “modes”  which  have  zero  displacement  at  the  upstream 
end. 


2.3  NONHOMOGENEOUS  TERM 


In  a  very  general  sense  the  ocean  space-time  field  may  be  thought  of  as  an  ad¬ 
dition  of  plane  waves  and  may  be  written  as 

A  V  *  i(w,  t  -  Tc*"x)  . 
v  -  2-  A^(t)e  k  (8) 

f 

where  k  is  a  two-dimensional  wavenumber  and  the  summation  is  overall  values  of  k. 

Because  we  are  primarily  concerned  with  towed  sensors,  we  must  account  for 
Doppler  effects..  If  we  are  towing  in  a  direction  parallel  to  the  negative  x  direction  at 
the  speed  U,  and  we  are  located  along  the  cylinder  at  the  position  x,  the  argument  of 
the  exponent  of  equation  (8)  becomes 

k  sin  ^Ut  -  k  sin  <^x  +  ui^t  ,  ^ 


where  lie  I  =  lie  I  and  $k  is  the  angle  of  k  with  respect  to  the  y  axis.  The  oceanic 
process  will  have  a  phase  velocity,  Cfc,  defined  from  the  dispersion  relation,  kCfc  = 
a>k-  Equation  (8)  thus  becomes 


V(x,t)  -£A£<t)eil(ulk  +  kxU)t  -  kxxJ  * 

* 


(10) 
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where  c* >\r  +  kxU  is  a  Doppler-shifted  angular  frequency  and  kx  is  the  x  component 
of  k.  Thus,  in  the  translating  array  coordinate  system  the  ocean  currents  appear  like 
plane  waves  traveling  in  the  positive  x  direction  and  having  an  alongtrack  wave 
number,  kx,  and  a  frequency,  kx  (Cfc  +  U).  Typically,  Cjc«U.  These  can  be 
superimposed  in  the  usual  Fourier  sense. 

We  next  drop  the  x  subscript  on  the  wavenumber  and  recognize  that  the  wave- 
number,  k,  is  the  alongtrack  wavenumber.  The  argument  of  the  exponent  in 
equation  (10)  may  be  written  as 

+  -  “-[(t  +  2)  T  -  !  ] 

which  is 

k[(Ck  +  1)  T  -  51  . 

where  the  “hat”  denotes  a  nondimensional  variable.  Henceforth,  in  this  report  the 
“hat”  will  be  dropped  and  k,  Cfc,  and  a>  will  be  nondimensional. 


3.  APPROXIMATE  SOLUTIONS 


Approximate  solutions  that  apply  to  low-  and  high-frequency  regions  are  ex¬ 
plored  in  this  section  because  they  lead  to  a  physical  understanding  of  the  process. 


3 . 1  LOW-FREQUENCY  APPROXIMATION 

Using  a  perturbation  technique  to  solve  equation  (1),  it  was  shown2  that  an 
important  solution  to  the  forced  vibration  problem  occurs  when  the  forced 
vibration  frequency  is  small, 3  that  is,  when  wL/U  <  1.  A  low-frequency  vibration 
convects  back  along  the  cylinder  relatively  unattenuated.  At  higher  frequencies 
significant  viscous  damping  and  structural  stiffness  take  place.  Thus,  in  actual 
towed  systems  typically  subjected  to  a  range  of  excitation  frequencies,  the  principal 
dynamics  take  place  at  the  lowest  frequencies.  A  similar  result  would  be  anticipated 
in  our  present  non-homogeneous  case,  i.e.,  equation  (7). 


Reference  2  shows  that  in  the  low-frequency  approximation  the  inertial  force, 

82n  82n 

312  +  28  • 

and  the  second-order  part  of  the  tension  stiffness,  (a  +  b)  d2rj/342,  are  negligible 
relative  to  the  hydrodynamic  forces  and  the  first-order  part  of  the  tension  stiffness, 
i.e., 

cln 

C3C  * 

Note  that  the  two  parts  of  the  tension  force  are  a  consequence  of  the  linearly 
decreasing  tension  of  a  towed  sensor.  We  see  that  in  dimensional  quantities 
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T(x) 


a2y  .  3T(x) 
3x2  8x 


(ID 


where  dT(x)/dx  is  a  constant  and  T  is  tension.  Thus,  the  low-frequency  ap¬ 
proximation  to  equation  (7)  is 


|2  +  c  !f-v 

3t  o  35 


(12) 


where 


c 

o 


References  6  and  9  show  experimentally  that  rtct/cn  is  of  order  10"! . 


Following  reference  7,  we  let  u(|,T,rj)  =  constant  be  an  integral  of  equation  (12). 
Then 


3u 

35 


3u 

3n  35 


V  u  V  H 


and 


Substitution  of  equation  (13)  into  (12)  yields 

+  c  +  V|£ 
3t  o35  3n 

which  may  be  written  as  a  vector  product 


3u  3u  3n 

3t  +  3n  3t  u 

(13) 

o  . 

(14) 

•  o  , 

(15) 

where  i,  j,  and  k  are  unit  base  vectors  in  ^Cartesian  coordinate  system  and  A  is  the 
gradient  operator.  Thus  the  vector?  +  c^  -  Vk  is  perpendicular  to  the  normal  to  u 
and  is  thus  tangent  to  a  solution  surface. 


The  vector  indicates  the  direction  in  l.r.rj  space  taken  by  the  solution.  The  path 
taken  by  a  “particle”  in  the  solution  space  is  called  the  characteristic  of  the 
equation.  The  equation  for  the  characteristic  is 

dT  .  f  -  f  .  (16) 

o 

which  defines  two  ordinary  differential  equations  that  have  the  integrals 

^(5. T,n)  -  and  u2(5,T,n)  -  c2  . 

Any  intersection  of  these  two  surfaces  is  the  characteristic  curve.  The  intersection 
occurs  1®  whenever 


u!  -  f(u2)  , 

where  f( )  is  an  arbitrary  function.  The  integrals  of  c0dT  =  d£  and  -Vd4  =  c0drj  are 

•  5  -  cqt  -  cl  (18a) 
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and 


u2  ■  n  + 


f  J 


d?V 


(18b) 


The  intersection  of  U}  and  U2  surfaces  gives  the  solution 


-ff 


d£V  -  f(g  -  CqT)  , 


(19) 


where  f( )  may  be  evaluated  from  the  upstream  boundary  condition,  rj(0,T)  =»  g(r). 

In  the  previous  section  we  saw  that  V  may  be  written  as  the  sum  of  plane  waves. 
Thus  we  are  interested  in  the  integral 


f 


d?V  "  "ik"  eikt(Ck  +  l)  T  "  ^  +  K(T)  for  k  <  1  ,  (20) 


where  k(T)  is  the  integration  constant.  The  solution  we  seek  then  becomes 

ik[(c.  +  1)t  -  a 


n(?,T)  -  f(?  -  c  t)  +C*' 


ik? 


-  K(T) 


where  $  =  |v|/cQ  . 

In  order  to  meet  the  upstream  boundary  condition,  that  is, 

n(0,T)  -  g(T)  , 


we  set 


and 


f(  )  -  g(  ) 
+ik(ck  +  1)  t 


K(t)  -  - 


ik 


<*  • 


With  these  conditions  the  solution  becomes 


i[k(ck  +  1 ) T  -  kCl 

n(?,t>  -  8(5  -  coT)  +  ?<J>  ^ ^ -  (eik^  - 


1) 


for  k<  1 .  For  k  <  1  it  is  useful  to  make  the  approximation 

-ik? 


±  ~ 

ik?  \ 


ik? 


Thus  to  first  order,  the  solution  is 


ik?  -  + 


-) 


H(?.T)  -  g(?  -  c  t)  +  ?$ei((u)k  +  U)T  ~  for  k  <  1  . 


(21) 
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The  solution  contains  two  parts:  (1)  the  boundary  condition  propagating 
downstream  nondispersively  with  the  velocity  c0,  and  (2)  the  nonhomogeneous 
forcing  functions  resulting  in  an  additive  solution  that  is  proportional  to  the 
magnitude  of  the  angle  of  the  local  fluid  velocity,  |v|/c0.  Physically,  this  is  a 
consequence  of  the  cylinder’s  natural  tendency  to  align  itself  with  the  flow,  which  in 
this  case  is  the  orthogonal  addition  of  U  and  $.  Actually  the  magnitude  of  the 
cylinder  slope  with  respect  to  the  x  axis  is  not  proportional  to  $/U  as  might  be 
thought,  but  rather  is 

a  A 

sin  k  [<ck  +  l)Ut  -  xl  .  (22) 

o 

Physically,  the  actual  angle  is  decreased  by  the  impedance  of  the  slope-dependent 
tension  stiffness  (equation  (1 1))  term,  cQ.  Because  the  cylinder  nearly  aligns  with  the 
flow,  the  displacements  are  the  spatial  integral  of  the  slope,  which  enters  the 
wavenumber  in  the  denominator. 

The  f(|  -  cqt)  term  in  equations  (19)  and  (21)  is  the  “water  pulley”  effect 
described  by  Kennedy.^ 

Noting  from  equations  (16)  and  (22)  that  the  space-time  slope  of  the  cylinder  is  - 
V/cQ  we  may  rewrite  equation  (19)  as 


n(£.T)  -  g(£  -  cqt)  -  J  d£8(£,t)  ,  (23) 

where  we  use  a  small  angle  approximation  for  6  which  is  explicit  in  the  derivation  of 
equation  (1).  Thus  we  see  that  at  low  frequencies  the  particular  integral  of  equation 
(7)  results  from  the  cylinder  locally  aligning  with  the  flow  so  as  to  minimize  the 
hydrodynamic  resistance.  The  shape  of  the  system  is  then  the  spatial  integral  of  the 
slope  by  definition. 

3.2  HIGH-FREQUENCY  APPROXIMATION 


Returning  to  the  PDE  of  equation  (7),  we  again  summarize  the  coefficients  as 
a  -  8  -  28cte  , 


b  -  2Bcte  , 

c  ■  2e8  (e  +  c  /it)  ,  and 
c  n 

c 

d  -  28e—  . 

IT 


(24) 


A  significant  simplification  results  from  exploiting  the  fact  that  the  cylinders  of 
interest  have  quite  high  aspect  ratios,  that  is,  e  is  large.  A  further  simplification 
results  from  neglecting  ctn/cn  relative  to  unity.  As  noted  previously,  ctn/cn  is  of 
order  10*  1 .  Thus,  if  t  -*  »  equation  (7)  becomes 


(1  -  5) 

ar 


Cn  3n  _  cn  3ri 
TTC  3£  ~  TTC  3t 

W  U 


(25) 
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From  the  discussion  following  equation  (10)  we  see  that  high  frequencies  are 
equivalently  high  wavenumbers  because  u>  at  kU.  Thus,  a  high-frequency  ap¬ 
proximation  is  equivalently  a  high  wavenumber  approximation.  Again  following 
the  procedures  of  reference  2,  we  make  a  WKB-like  approximation  in  assuming  that 
d^rj/dl2  changes  rapidly  during  small  changes  in  (1  - 1).  If,  therefore,  (1  - 1)  is 
treated  as  a  constant,  then  equation  (25)  becomes  a  linear  PDE  with  constant 
coefficients.*  Because,  in  this  approximation,  the  cylinder  takes  on  numerous  cycles 
and  our  motivation  is  only  to  identify  physical  mechanisms,  we  simplify  the 
problem  by  neglecting  the  boundary  conditions  and  making  the  cylinder  infinitely 
long.  If  we  let 

n(S,T)  -  r(T)e"ikC  ,  (26) 

then  equation  (25)  becomes 


(1  -  «  +  ik  if 


c  T 
n 

Tct  Co 


r(l) 


Cn  dr(r) 
irct  di 


(27) 


The  Laplace  transform  of  equation  (27)  for  the  case  of  zero  initial  conditions  is 

R(s)  (s  +  Q)  -  -V(£,s)  eikC  (28) 


where  R  and  V  are  transforms  of  r  and  V,  s  is  the  transform  variable,  and 

q  -  -ikcQ  +  k2  ^  d  -  C)  .  (29) 

n 

Directing  our  interest  to  the  plane  wave  case,  that  is,  V(|,t)  *  V  exp  i  (tor  -  k£), 
we  have 


R(s) 


-  V 


(s  +  Q)  (s  -  iu>) 


The  inverse  transform  is 


r(T)  -  — | V I  A(e~^T  -  eiwT)  , 


where 


or 


A  -  (Q  -  iu>) 


-1 


TTC 


-1  ^ct  2 

A  — -  (1  -  5)k  +  i(kc  -  u) 

c  o 

n 

(kco  •  w)  WET 
1  +  i  - - - 5— s- 

(1  -  £)kz 


-  (1  -  C)k2 


(30) 


(31) 


(32) 


*It  was  found  in  reference  3  that  this  approximation  leads  to  an  adequate  model 
for  k  >50  and  4  <  0.98. 
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For  our  current  large  k  case. 


(kco  -  M>  157 
(1  -  Ok2 


<  < 


(33) 


Therefore, 


_ n_ 


A  — 


k  (1  -  C) 


1  -  i 


^7  (kco  - 
(1  -  C)k2 


c 

_ n_ 

TTC_ 


k  (1  -  O 


where 


Can  ip 


rr-  (kc  -  til) 

7TCt  O 

k2 ( 1  -  O 


Finally  we  write 


n(C,x) 


_ n 

TTC 


k  (1  -  5) 


t_  ei4» 


i(u)T  -  kC) 


-k2  c  (1  "  C)T  i(kc  T  -  kO 
—  e  f  n  e  o 


(34) 


(35) 


(36) 


for  k  >50  and  $  <0.9. 

Because  of  the  large  values  of  k  necessary  for  this  approximation,  to  hold  the 
second  term  in  the  brackets  rapidly  becomes  negligible,  and  the  solution  is  a  wave 
propagating  downstream  at  the  tow  speed. 

The  amplitude  of  the  wave  is  exactly  the  ratio  of  the  normal  hydrodynamic  force 
applied  by  the  crosscurrent  on  the  cylinder  to  the  tension  stiffness  of  the  cylinder 
tensioned  by  tangential  hydrodynamic  drag  for  a  sinusoidally  perturbed  cylinder. 
Using  a  hat  to  denote  dimensional  quantities  we  show  the  above  statement  by  noting 
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that  the  normal  force  to  the  cylinder*  is  l/2<?DcnUV,  and  the  tension  stiffness  is 
T(x)d2y/  3x2  where  T(x)  »  l/2eDctLlj2(l  -  x/L)  .  Thus,  if  the  induced  shape  is 
assumed  sinusoidal  with  maximum  amplitude  A  and  wavenumber  k,  then  the  two 
forces  will  balance  when 


which  in  nondimensional  form  is  exactly  the  amplitude  of  equation  (36).  Thus  the 
amplitude  of  the  wave  is  determined  by  string-like  forces,  but  not  the  propagation 
velocity.  The  significant  amplitude  attenuation  caused  by  tension  stiffness  is 
reduced  as  (1  - 1)'1  reflecting  the  tension  reduction  as  one  approaches  the  free  end. 


4.  COMPLETE  SOLUTION 

The  discussion  in  this  section  is  limited  to  the  solution  to  equation  (7c)  because 
the  solution  to  equation  (7b)  has  been  detailed  in  reference  1.  The  subscript  (ni)  of 
equation  (7c)  will  be  dropped  to  increase  notational  efficiency. 

In  solving  equation  (7c),  we  make  the  fundamental  assumption  that  the  solution 
of  »t(4,t)  is  separable  in  space  and  time,  and  the  time  dependence  is  harmonic,  that 
is, 

n(?,x)  -  s(Oeia>T 
v(C.T)  - 


where  co  is  noted  to  be  nondimensional  and,  in  general,  complex,  and  i  denotes  the 
imaginary  number.  Note  that  after  the  substitution  the  primes  are  dropped. 

Substituting  the  above  equation  into  equation  (7c),  along  with  the  transformation 
of  the  independent  spatial  variable  (z2  =  a  +  b|),  yields  the  ordinary  differential 
equation  (ODE) 

^  +  b2s  -  c  v  for  0  <  z2  <  a  (38) 

.  2  z  dz  o  o  v  ' 

dz 


with  the  boundary  conditions 

2 

S(z)  ■  0  for  z  ■  a  , 

|s(z)|  <  00  for  0  <.  z2  <_  a  , 


•This  form  was  originally  used  by  Paidoussis  in  reference  4  and  more  recently  by 
Ketchman  in  reference  1 1 . 
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where 


£  (c  +  2iuB)  -  1  , 


bo  "  dd  “  »  ant^ 

b 


4d 
.  2 


Equation  (38)  is  recognized  as  a  modified  form  of  Bessel’s  equation.^ 
It  is  useful  to  make  the  transformation 

w(z)  -  S(z)zp  , 

where 

a  -1 
o 


(39) 


After  some  algebra  the  governing  ODE  becomes 

£w  +  X2zw  ■  cqzp+*V(z)  ,  (40) 

where  C  is  the  Sturm-Liouville  operator1® 

t.<U)  .£ 

*  dz  \  dz  /  z 

containing  the  coefficients  of  Bessel’s  equation  and  equals  bg  .  The  boundary 
conditions  for  this  problem  are 


v(z)  ■  0  for  z  -  a 


1/2 


(41) 


and  w(z)  is  finite. 

Much  has  been  written  of  the  Sturm-Liouville  problem.  An  important  aspect  of 
this  classical  ODE,  pertaining  to  our  present  problem,  is  that  because  z  vanishes  at 
the  upstream  end  of  the  interval  (i.e.,  z  =  0),  and  the  solution  vanishes  at  the  other 
end  of  the  interval  (i.e.,  w(al/2)  =  0),  then  the  Eigenfunctions  of 

£  w  +  A2zv  -  0  (42) 

n  n  n 

are  orthogonal  with  respect  to  the  weighting  function  z  over  the  interval  0  to  a  1/2, 
and  are  Eigenvalues  that  satisfy  the  boundary  conditions. 

Thus, 
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*1 

zw  w  dz  -  0  if  n  +  m  .  (43) 

n  m 

Jo 

Actually  this  differs  somewhat  from  our  problem  in  that  the  interval  is  a"^  + 
b/a)  to  a1^,  not  0  to  a1^.  This  is  a  consequence  of  the  switching  of  equation  (1) 
from  the  hyperbolic  PDE  to  an  elliptic  PDE  quite  near  the  downstream  end.1  All  of 
the  above  mathematics  apply  only  in  the  hyperbolic  region.  The  switchover  point1 
occurs  when 

Z  -  a/b  .  (44) 

Thus,  if  a/b  is  unity  the  switch  does  not  take  place  until  the  end  of  the  cylinder  is 
encountered  and  the  interval  of  z  becomes  0  to  1 .  Writing 

-a/b  -  1  -  l/2ct£  (45) 

we  see  that  this  condition  is  nearly  approached  for  long  cylinders.  We  assume  that 
this  condition  is  met. 

Because  equation  (42)  is  the  classical  Bessel  equation,  the  Eigenfunctions  are 
Bessel  functions  of  the  first  and  second  kind.  Bessel  functions  of  the  second  kind 
will  be  singular  in  the  interval  of  interest,  thus  they  are  excluded  by  the  finite 
condition1  of  equation  (41).  Therefore,  the  complete  solution  to  equation  (42)  with 
the  boundary  conditions  of  equation  (41)  is 

w(z)  -  £  bnJp<*nZ>  for  0  1  z  1  “1/2  »  (46) 


where 


P 


(47) 


and  knal/2  are  the  zeros  of  the  Bessel  function.  Following  reference  1  we  avoid  the 
complication  of  complex-order  Bessel  functions  by  again  requiring  a  long  cylinder, 
long  enough  so  that 

|w|/(cte)  <  <  cn/irct  • 


With  this  condition  met  all  An  are  real  numbers.  The  coefficients,  bn,  are  found 
from  initial  conditions. 


Because  of  the  orthogonality  of  the  Eigenfunctions  the  nonhomogeneous 
equation  (40)  can  be  readily  solved  by  assuming  a  solution  of  the  form 

v(z)  -  £  anwn  ,  (48) 

n 
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which  can  be  substituted  into  the  equation  to  be  solved  yielding 


E(£a  w  +  X^za  w  )  «*  c  z^^Vfz) 
v  n  n  n  n  o 


(49) 


Using  equation  (42)  we  simplify  the  above  equation  by  writing 

_  2  ,  c/Vz) 

£  (X2  -  X2)a  w  -  — — - - 

n  n  n  z 


(50) 


Next  we  express  the  right-hand  side  of  the  above  equation  as  a  series  of  orthogonal 
functions,  that  is, 


n  n  z 

n 


(51) 


We  evaluate  the  coefficients  by  multiplying  both  sides  of  equation  (51)  by  wm  and 
integrating  over  the  interval,  that  is, 

1/2  .1/2 


/•a 


n  J 


f  +1 

dzzA  w  w  *  c  I  dzz^+^V(z)w 
n  n  m  o  I 

0  JO 


m 


Because  of  the  orthogonality  relation  of  equation  (43)  we  may  write 


.1/2 


A  I  dzzv 
n 


1/2 

,  f 

'  -  C  1 

n  o  1 
JO 


dzz^+*V(z)w  , 
n 


where  the  integral  of  the  left-hand  side  has  been  evaluated10  as 

L1/2 

dzz[J  (X  z)]2  -  a/2[J  .(Xa1/2)]2 

p  n  p~*  n 


Substitution  of  equation  (51)  into  (50)  and  equation  (54)  into  (53)  yields 

1/2 


'o  dzzzPv(z)Jp(*nz) 

yJ2  (X  a1/2)(X2  -  X2) 
2  p+1  n  n 


X2  -  X  2 
n 


(52) 


(53) 


(54) 


(55) 
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Finally,  after  combining  equations  (55),  (48),  and  (54)  we  write  that 

1/2 


S(z) 


£ 

n 


2 /a 


T2  t\  *1/2^ 

Jp+l(Xna  5 


yypV(y)j  (X  y)dy 
P  n 


z  PJ  (X  z) 
p  n 


,  (56) 


where  we  note  that  the  integral  is  the  Fourier-Bessel  Series  of  zPV(z)  and  z‘PJp(Anz) 
are  the  normal  modes.  This  equation  can  be  written  more  compactly  as 


V(y)J  (X  y) J  (X  z) 
_ p  n  p  n 

,,  1/2.  p.,,2  .2. 

(X  a  )z  (X  -  X  ) 
n  n 


dy 


(57) 


where  the  problem’s  Green’s  function  or  impulsive  response  is 


G(z,y) 


£ 


(y/z)pyJ  (X  y)J  (X  z) 

_ E _ 2 _ P  n 

a  f2  1/2. ^..2  .2. 

tjp+i(V  )(X  -  V 


(58) 


Plots  of  the  normal  modes,  z*Pjp(Anz),  normalized  to  unity  at  4  =  1  are  shown 
versus  cylinder  position  in  Figure  3.  Elements  of  the  approximate  solution  of 
equations  (21)  and  (36)  are  evident  in  the  plot.  The  lower-order  modes  show  the 
approximately  linear  growth  of  equation  (21)  and  the  higher-order  modes  show  the 
inverse  tension  (i.e.,  (1  - 1)'1)  behavior  of  equation  (36).  As  expected,  the  complete 
solution  integrates  the  two  physical  characteristics  that  are  highlighted  in  the  ap¬ 
proximate  solutions,  that  is,  (1)  the  cylinder  attempts  to  align  with  the  flow  leading 
to  a  linear  growth  in  the  transverse  displacement,  but  (2)  the  cylinder  is  impeded  by 
cylinder  tension  stiffness  that  increases  as  one  moves  forward  along  the  cylinder  and 
increases  with  wavenumber  increase.  The  frequency  or  wavenumber  dependence 
will  be  determined  by  the  magnitude  of  the  various  modes  which,  in  turn,  are  ob¬ 
tained  from  the  Fourier-Bessel  spectrum  of  zPV(z).  Unfortunately,  the  spectrum  of 
a  planewave  is  rather  broad,  so  the  series  of  equation  (56)  does  not  rapidly  con¬ 
verge.  In  fact,  equation  (56)  is  not  a  numerically  efficient  approach  to  obtaining 
solutions.  A  method  for  generating  numerical  solutions  is  discussed  in  the  next 
section. 


The  wavenumber,  or  equivalently  frequency,  dependence  of  the  system  is 
analytically  determined  by  the  wavenumber  dependence  of  the  Fourier-Bessel 
spectra  discussed  above  and  the  frequency  dependence  of  the  (A^-A^)'1  term  in 
equation  (56).  Using  the  parameters  defined  in  equation  (38)  we  write 


X 


2 


n 


where  dn  is  the  n1*1  real  zeros  of  the  p^  order  Bessel  function.  Employing  the 
assumptions  discussed  following  equations  (45)  and  (47)  and  the  parameters  defined 
in  equation  (38)  we  write 
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•'S' 


'?!* 

$ 

$ 


X2  -  X2  .  -14  ^ 
n  ire 


2Bcte 


j  4c  /  ir62  \ 

-  X2  s - 2-  [ioj  +  - ~ 

n  »ct  V  86cnc  /  ' 


Therefore,  we  conclude  that  (A^  ‘  ^n2)'1  is  a  simple  pole  reflecting  viscous 
damping  of  transverse  motion.  In  other  words,  the  “velocity”  or  “dashpot” 
damping  of  equation  (3)  leads  to  a  (icu)*1  dependence. 


5.  NUMERICAL  SOLUTIONS 

The  analytical  solution  of  section  4  has  given  a  better  physical  understanding  of 
the  process  by  identifying  the  constitute  modes  of  behavior  of  the  system.  The 
method  is  not,  however,  a  numerically  efficient  way  of  realizing  specific  space-time 
array  shapes.  In  this  section  an  efficient  method  of  obtaining  a  numerical  solution  is 
described. 

A  solution  to  the  forced  vibration  problem  of  equation  (7b)  using  a  finite  dif¬ 
ference  scheme  in  conjunction  with  powerful  numerical  ordinary  differential 
equation  methods  is  described  in  reference  8.  The  method  of  reference  8  has  been 
expanded  to  solve  the  current  problem  of  interest  summarized  in  equation  (7c). 

Reference  8  details  the  procedures  which  are  merely  summarized  here  for  reader 
convenience. 

We  proceed  by  applying  the  second-order  central  and  backward  discre  ization  to 
the  partial  differential  equation  (7)  in  the  |  direction.  This  yields  the  form 


+  (a  +  b£ 


b{il(VrVV.)  +c  , 
f  £(n»  -  Vi)  +  d  IT  ■  -« 


+  h  <n«  -  Vi> 


where 


5_  =  mh,n  =  and  V  =  V(E  ,x) 

m  m  m  m  ra 


Equation  (59)  is  a  difference  equation  representing  a  system  of  second-order 
ordinary  differential  equation  and  is  an  approximation  of  equation  (7).  This  is 
converted  to  a  system  of  first-order  ordinary  differential  equations  by  substitution 
of 


into  equation  (59).  After  some  manipulation  the  first-order  system  may  be  written 
in  matrix  form  as 


£  N  -  A(C)n  +  g(C,T,n)  (61) 

where 


W1 

^  w  +  £  n  -dV 
h  o  h  o  1 

N  -  , 

n2 

and 

g  " 

0 

W2 

-dv2 

Vi 

WM-1 

0 

a  +  bt,  i 

h2  -%-! 

and  M  is  the  number  of  spatial  points. 


The  g  vector  contains  only  the  boundary  conditions  and  nonhomogeneous  for¬ 
cing  particular  to  the  present  problem.  The  A(4)  matrix  contains  no  forcing  terms  at 
all  and  thus  represents  a  known  linear  transformation.  The  difference  between  the 
above  equation  and  equation  (3.10)  of  reference  8  is  the  inclusion  of  the 
nonhomogeneous  terms.  These  terms  are  inputs  to  the  problem  in  the  same  sense  as 
boundary  conditions  are  inputs.  Thus,  the  solution  of  equation  set  (61)  can  proceed 
exactly  as  it  did  in  reference  8  using  the  generalized  Adams-Bashforth  method  —  a 
highly  efficient  procedure  for  this  type  of  problem. 

Application  of  the  upstream  (or  top)  boundary  condition  of  equation  (7c)  makes 
wo  =  *7o  =  0  in  ^e  g  vector.  However  the  boundedness  condition  of  equation  (7c) 
does  not  adequately  specify  r)M  which  is  an  explicit  value  of  the  downstream  (or 
bottom)  boundary  condition.  This  is  an  important  difficulty  associated  with  the 
application  of  these  techniques  to  this  particular  problem  (refer  to  the  discussion  of 
boundary  conditions  in  section  2  of  this  report).  However,  because  the  existence  of 
the  nonhomogeneous  term  in  the  present  problem  in  no  way  modifies  the  down¬ 
stream  boundary  condition,  we  may  use  the  approximate  ad  hoc  condition  derived 
in  reference  8.  The  approach  was  to  integrate  the  transverse  momentum  equation 
(equation  (1))  over  a  short  tapered  end  which  is  attached  to  the  free  end  in  order  to 
generate  the  required  boundary  condition.  The  result  of  this  approach  is  the  con¬ 
dition 


In  a.  in 

3T 


0 


(62) 


which  is  shown  in  reference  8  to  adequately  solve  the  problem.  Pnysically  this 
condition  represents  a  "radiation  condition”  which  reflects  no  energy. 
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Figure  4a.  Spatial  Cylinder  Shapes 


Figure  4.  Cylinder  Shapes 
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As  was  the  case  in  the  previous  two  sections,  we  are  interested  in  solutions  to  the 
“planewave”  case  discussed  in  section  2.  This  was  accomplished  by  numerically 
solving  an  initial  value  problem  having  the  following  conditions 


n(S.o)  -  0  for  0<  5  <  1, 

dnd(^,0)  -  0  for  0  <  1,  and 

V(£,t)  -  ei(“T  "  k5)  for  0  <  £  <  1  . 


(63) 


The  results  of  the  technique  are  displayed  in  two  ways:  (1)  transverse 
displacement  (rj)  versus  longitudinal  position  (|)  with  time  (t)  as  a  parameter,  and 
(2)  transverse  displacement  versus  time  with  longitudinal  position  as  a  parameter.  In 
the  first  type  of  display  (figure  4a),  each  sequential  time  increment  is  offset 
downward  in  a  “waterfall"  type  display.  These  displays  show  qualitatively  the  two 
characteristics  of  the  solution  discussed  in  the  previous  two  sections,  that  is,  (1)  the 
solution  appears  to  propagate  aft  with  unity  nondimensional  phase  velocity  which 
in  dimensional  quantities  is  the  axial  flow  velocity,  U;  and  (2)  the  envelope  of  the 
solutions  grows  in  the  downstream  direction.  The  latter  phenomenon  is  better 
illustrated  in  figure  4b  which  shows  that  the  time  periodic  steady  state  response 
increases  downstream,  that  is,  increasing  |. 


6.  SUMMARY 

It  has  been  shown  that  the  transverse  dynamics  of  a  long  flexible  cylinder  im¬ 
mersed  in  axial  fluid  flow  and  perturbed  by  a  space-time  nonaxial  flow  component 
can  be  modeled  by  a  linear  nonhomogeneous  PDE.  The  subject  momentum 
equation  is  related  to  a  previously  studied  forced  vibration  problem  that  was 
modeled  by  a  similar  homogeneous  PDE.  The  intent  of  this  work  is  to  model  the 
transverse  dynamics  of  a  cable/array  system  driven  by  crosscurrents,  and  to  relate  it 
to  a  previous  study  of  a  cable/array  system  excited  by  the  random  crosstrack 
meandering  of  the  towing  vessel. 

Three  types  of  analytical  solutions  were  explored  in  order  to  develop  a  physical 
and  mathematical  understanding  of  the  ocean-induced  transverse  motion  of  towed 
arrays.  The  following  physical  picture  emerges  from  this  effort: 

1 .  In  cable/ array  systems  having  aspect  ratios  (length  to  diameter)  greater  than 
lO^,  inertial  forces  can  be  neglected. 

2.  If  the  cable/array  system  length  is  much  shorter  than  the  alongtrack 
wavelength  of  the  oceanic  current  system,  then  hydrodynamic  forces  dominate 
structural  forces  and  the  cable/ array  locally  aligns  with  the  flow.  Transverse 
displacements  are  then  the  spatial  integral  of  the  cable/ array  angle  with  the  fluid. 
The  displacement  response  of  the  system  will  be  proportional  to  V/Uk,  and  these 
displacements  will  appear  to  propagate  aft  at  nearly  the  tow  speed. 
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3.  If  the  cable/array  system  length  is  much  longer  than  the  alongtrack  wavelength 
of  the  oceanic  current  system,  then  the  tension  stiffness  of  the  cylinder  is  the 
principal  impedance  to  the  transverse  ocean  forces  applied  to  the  cable.  The  con¬ 
sequence  of  this  fact  is  that  the  attenuation  of  the  transverse  dynamics  by  tension 
stiffness  is  very  much  a  function  of  the  spatial  position  on  the  cylinder  and  the 
alongtrack  wavenumber  of  the  ocean  current.  Because  the  tension  goes  to  zero  at  the 
free  end,  the  dynamics  near  the  free  end  are  larger  relative  to  upstream  positions. 
Displacement  response  of  the  system  will  be  proportional  to  (^/UKCn/nc^/k^L  - 
x),  where  cn/nct  is  a  measure  of  the  normal  to  transverse  hydrodynamic  drag.  These 
displacements  will  appear  to  propagate  aft  at  nearly  the  tow  speed  and  grow 
significantly  near  the  free  end. 

4.  Normal  hydrodynamic  drag  acts  as  a  “velocity”  or  “dashpot”  damper 
yielding  a  frequency-dependent  dissipation  term  proportional  to  (cu)'1. 

5.  The  magnitude  of  the  “tail”  dynamics  relative  to  upstream  positions  is  a 
function  of  wavenumber,  or  equivalently  apparent  frequency,  of  the  oceanic 
current  process.  At  small  wavenumbers,  where  hydrodynamic  forces  dominate 
structural  stiffness,  little  “tail  wagging”  will  be  observed.  At  large  wavenumbers 
the  envelope  of  the  “tail  wagging”  will  be  proportional  to  (L  -  x)*l . 

6.  The  actual  response  of  a  cable/array  system  to  ocean-forcing  is  determined  by 
the  Fourier-Bessel  spectra  of  a  function  of  the  spatial  distribution  of  the  crosstrack 
ocean-current  structure.  Such  spectra  need  to  be  studied  to  determine  the  actual 
spectra  of  ocean-induced  cable/towed  array  transverse  dynamics. 

7.  The  transverse  shape  of  a  cable/array  system  can  be  generally  expressed  as  a 
sum  of  “modes”  of  the  form  z-pJp(Anz).  The  convergence  of  this  sum,  for  the  case 
of  planewaves  and  impulsive  spatial  structures,  is  rather  slow,  so  this  procedure  is 
not  a  very  numerically  efficient  method  of  solution. 

8.  A  numerical  solution  to  the  subject  governing  equation  which  uses  a  finite 
difference  scheme  (i.e.,  “method  of  lines”)  in  conjunction  with  powerful  numerical 
ordinary  differential  equation  methods  (i.e.,  “generalized  Adams-Bashforth” 
method)  has  been  programmed  and  examples  run. 

There  are  three  general  areas  of  future  work  necessary  to  understand  the  ocean- 
induced  towed  array  dynamics  problem:  (1)  further  quantitative  comparison  bet¬ 
ween  analytical  and  numerical  procedures  is  necessary  to  give  credibility  to  the 
numerical  technique;  (2)  the  Fourier-Bessel  spectra  of  typical  high  wavenumber 
ocean  processes,  such  as  internal  waves  and  upper  ocean  turbulence,  are  needed  to 
understand  what  structural  modes  will  be  excited;  and  (3)  a  quantitative  comparison 
between  theory  and  existing  data  is  needed  to  give  credibility  to  this  modeling 
procedure. 
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